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We show that the low-energy sector of the S=l/2, antiferromagnetic Heisenberg model on a one- 
dimensional lattice of coupled tetrahedra consists of 2^ replica of the spectrum of the dimerized 
Heisenberg chain, where N is the number of tetrahedra. This provides a proof of the following 
properties: i) there is a residual ground-state entropy per spin equal to 2 1//4 ; ii) there is a singlet- 
triplet gap as long as the coupling between the tetrahedra is smaller than the internal one. These 
properties are compared to available results on the pyrochlore lattice. 
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It is by now quite clear that the low-energy properties 
of frustrated magnets can be very different from those 
of ordinary magnets with long-range order. In particu- 
lar, after the pioneering work of Majumdar and Ghosh 
on the zigzag chain |l]], it has been shown that several 
systems have a singlet-triplet gap in the magnetic spec- 
trum. This property is actually not specific of frustrated 
systems since it is also shared by ladders with an even 
number of legs || . There is an increasing evidence how- 
ever that frustration can have more specific consequences, 
like for instance low-lying singlets in the singlet-triplet 
gap. The first example was again the zigzag chain which, 
contrary to ladders, has a two- fold degenerate singlet 
ground-state |l| . But the spectrum can be more complex. 
For instance, there is a clear numerical evidence that the 
Heisenberg antiferromagnet on the Kagome lattice has an 
exponential number of low lying singlets below the first 
triplet excitation ||,|| . If analytical arguments have been 
put forward to explain this property in terms of coupled 
triangles ||, an analytical proof that this is the case is 
still lacking. Such properties are actually not limited to 
low-dimensional systems. In particular, the Heisenberg 
model on the pyrochlore lattice is a well known exam- 
ple of 3D frustrated model without long-range order at 
the classical level ||, and the S=l/2, quantum version 
is believed to have a very short correlation length and a 
singlet-triplet gap in the ground-state 0. There is also 
some indication that there are low-lying singlet states in 
the singlet-triplet gap 0] , but again there is no proof that 
this is indeed the case. 

In this paper, we show that these properties actually 
occur in a one-dimensional analog of the pyrochlore lat- 
tice, namely a one-dimensional array of coupled tetrahe- 
dra. The discussion can actually be carried out for the 
slightly more general situation of a one dimensional sys- 
tem of alternating spins and triangles (see Fig. defined 
by the Hamiltonian: 
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where S n ^ are spin 1/2 operators. Similar models involv- 
ing more than one spin at a given site have been studied 
^|-|lT| , but to our knowledge none of them dealt with tri- 
angles, and the properties that are described below are 
very specific to that case. 
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FIG. 1. System of alternating spins and triangles 

Ti = {(2, i), (3, i), (4, £)}. The notations for the couplings are: 
Ji for bonds inside triangles and J2 (respectively J3) between 
a triangle Ti and its left (right) spin neighbor Sij (Si,i+i). 

The most important step toward a solution of this 
Hamiltonian is to realize that the spins of the triangle of 
any unit cell i enter the Hamiltonian only through their 
sum Ti — S2,i + £3, i + £4,3 since the first term can be 
rewritten ( Ji/2) — 9/4). This observation has two 

consequences. First of all, it means that the total spin of 
the triangle of any unit cell, which can take the degen- 
erate value 3/2 or the twofold degenerate value 1/2, is a 
good quantum number. Secondly, it shows that the ad- 
ditional quantum number needed to specify the states in 
the subspace of total spin 1/2 - for instance the chirality 
- does not enter the Hamiltonian. 

So the eigenvalues of the problem of Eq. (|l|) are the 
same as those of the following family of effective Hamil- 
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tonians H({Ti}): 

= ^(J 2 S^.fi + J 3 fi.S lji+1 ) 

i 

+ (J 1 /2)^(i; 2 -9/4) (2) 

i 

where 7^ 2 = T^T, + 1) and T t = 1/2 or 3/2, the de- 
generacy of the eigenvalues of H({Ti}) for the original 
problem being equal to 2™ 1 / 2 , where n 1 / 2 is the number 
of triangles with total spin 1/2. 

The twofold degeneracy associated with each triangle 
in a doublet state can actually be understood in a more 
direct way. Let us consider the eigenstates of the problem 
obtained by considering only one spin in the triangle of 
unit cell i, say 5*2, i- Then the wave functions obtained as 
the product of the eigenstates of that new problem with 
the singlet constructed out of £3^ and £4,1 are eigenstates 
of the original problem since all the extra couplings con- 
nect a single spin to both ends of that singlet. Besides, 
£3^ and £4^ forming a singlet, these eigenstates corre- 
spond to Ti = 1/2. Finally, the wave functions obtained 
by putting the singlet on the three possible bonds of one 
triangle are not linearly independent but generate a space 
of dimension 2, as for a single triangle. 

So the problem has now been split into different sectors 
corresponding to the values of the TVs. For clarity, we 
will continue the discussion in two specific situations: i) 
Alternating spins 1/2 and triangles, where the analysis 
is particularly straightforward; ii) Coupled tetrahedra, 
which is physically more relevant as an analog of the 
pyrochlore antifcrromagnet. 

I. ALTERNATING SPINS 1/2 AND TRIANGLES 

This case corresponds to J 2 = J3. For clarity, we re- 
name the parameters J\ = J and 3 2 = J3 = J'- So 
we are dealing in this section with triangles of strength J 
coupled to spins 1/2 with 3 exchange integrals of strength 
J'. In the limit where J' / J goes to zero, the spectrum 
of each sector Ti is completely degenerate, and the en- 
ergy is given by the sum of the energies of the triangles 
E({Ti}) = -(3/2)(N - J2i T i), wh ere N is the number 
of unit cells. In that limit, the lowest energy states are 
those obtained when all the triangles are in a doublet. 
Perturbative arguments show that this will remain true 
up to a certain value of the ratio J'/ J. To study that 
problem quantitatively, we have performed exact diago- 
nalizations of finite clusters for the Hamiltonians of Eq. 
(||). The results are given in Fig. ||. 

Finite size analysis. It is well known that the ground 
state energy per site of the spin 1/2 Heisenberg model 
scales like cr, = — A/L 2 where = 1/4 — In 2. So 
the ground state energy per site in units of J for TL{{Ti = 
1/2}), which is nothing but the standard Heisenberg 
model with coupling J' up to a constant of order J, scales 



like ef> = -3/16 + ( J'/2J) eoo - A/L 2 where L = AN is 
the total number of sites of the system. Let us focus now 
on e ( l\ the GS per site of H({T^ l0 = l/2,T io = 3/2}), 
and let us denote by the thermodynamic limit of the 
energy difference between the GS energy of the spin 1/2 
Heisenberg model and the GS energy of the same model 
with one spin 1/2 replaced by a spin 3/2. Since there 

are no 1/L corrections for e£ , one should expect the 
following scaling for e^ 1 -*, 

ei^^+A^/L + Otl/L 2 ) (3) 

where Aoo = (3 — (J'/J)S 00 )/2 is the thermody- 
namic limit of the energy difference between the GS of 
H({T^ io = 1/2, T l0 = 3/2}) and ft({T 4 = 1/2}). Nu- 
merical simulations up to L — 28 with one spin 3/2 show 
that the scaling (|J) is very well verified and allowed us to 
extract the gap (see Fig. |^) which, as expected, is linear 
in J', and 5^ ~ 5/2. 
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FIG. 2. Mam figure. Ground state energy per site of 
H({Ti^ = 1/2, T = 3/2}) for a L = 16 cluster (i.e. 4 tetra- 
hedra) for different values of the number of triangles Ti = 3/2: 
(•) n-3/2 = 0, (■) n 3/2 = 1, (A) n 3/2 = 2 neighboring trian- 
gles, (♦) n s /2 = 2 next-neighboring triangles, {■<) n 3 / 2 = 3, 
(▼) n 3/ / 2 =4. Inset. Thermodynamic limit of the energy dif- 
ference between the GS of H({T^ io = 1/2, T io = 3/2}) and 
Tt({T = 1/2}). Up to J' c ~ 1.2J the low energy physics of 
the model is given by the spin 1/2 Heisenberg model. 

Discussion. It is then clear that the ground state re- 
mains in the sector {Ti — 1/2} up to J' c ~ 1.2 J, and in 
particular for the isotropic point J' = J. When this is 
the case, the low energy physics is given by the Hamilto- 
man of Eq. @ with all T,'s equal 1/2, which, as we saw, 
is nothing but the one-dimensional, spin 1/2 Heisenberg 
model with coupling J'. The only difference is that we 
now deal with 2 N replica of the spectrum. All properties 
of interest can be deduced from this mapping. In partic- 
ular, there is no singlet-triplet gap in the spectrum, the 
elementary excitations are deconfined spinons, and there 
is a residual entropy per spin equal to 2 1 / 4 . 
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II. COUPLED TETRAHEDRA 

This case corresponds to J2 = J\. For clarity we use 
the notations J2 = J\ = J and J3 = J' in the follow- 
ing. So we are now dealing with tetrahedra of strength 
J coupled by bonds of strength J'. Let us start from 
the limit J' = 0. In that case, the groundstate is ob- 
tained by putting each tetrahedron in its groundstate. 
Since the ground state of a tetrahedron is twofold degen- 
erate with energy —(3/2) J, the groundstate is 2 Ar -fold 
degenerate with energy — (3/2)NJ. The low lying ex- 
cited states correspond to putting one tetrahedron in its 
first excited state with an energy cost of J. With re- 
spect to the general classification of the states proposed 
at the beginning of the paper, the situation is slightly 
more complicated than in the previous case. First of 
all, the groundstate manifold in the limit J' = 0, of 
dimension 2 , contains only some of the eigenstates of 
Ti({Ti ~ 1/2}), the Hilbert space of that Hamiltonian 
being of dimension 2 2N . Besides, the manifold of the 
first excitation of energy J contains states correspond- 
ing to different effective Hamiltonians. This comes from 
the fact that the triplet excitation of a tetrahedron is 
threefold degenerate, and that only two of them can be 
constructed with a given triangle being a doublet, the 
other one corresponding to the groundstate of a spin 1 /2 
coupled to the spin 3/2 state of a triangle. So it is no 
longer clear a priori that the low energy sector can be de- 
scribed by only one of the effective Hamiltonians of Eq. 
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FIG. 3. Extrapolated values of the GS energy difference 
between H{{T ilHo = 1/2, T io = 3/2}) and H({Ti = 1/2}) (.) 
and singlet-triplet gap of TL({Ti — 1/2}) (o). Up to large 
values of J' including J' — J, the effective hamiltonian of the 
system is H({Tj = 1/2}) and the low energy properties are 
those of the spin 1/2 dimerized Heisenberg model. 

To address this point, we have again resorted to ex- 
act diagonalizations of finite clusters, and we have de- 
termined which effective Hamiltonian gives the ground- 



state and the first triplet excited state as a function of 
J'. The results are given in Fig. ||. As in the previous 
case, the groundstate is given by H({Ti — 1/2}) regard- 
less of the size of the system up to very large values of 
J', and in particular beyond the isotropic limit J' = J. 
For the first excited state, the situation is slightly more 
involved since H({Ti — 1/2}) is now equivalent to the 
dimerized Heisenberg model and has a gap if J' =/= J 
p"2| . For small enough J 7 , the first excited state is also 
given by H({Ti = 1/2}) for all sizes. However, this is no 
longer true if J' is beyond a certain value J' C (N) which 
increases with the size of the system. Since we are inter- 
ested in the properties of the system in the thermody- 
namic limit, we have performed a finite size analysis of 
the first excited state of Ti^T, = 1/2}) and of the GS of 
H({T,^ l0 = 1/2, T io = 3/2}) (see results in Fig. |). 

Finite size analysis. Contrary to the situation dis- 
cussed in the previous section, there are exponential cor- 
rections to the finite size energies if J' ^ J because there 
is then gap in the spectrum of H({Ti = 1/2}). Namely, 
with the same notations as in section I, 

e^^e^+A^/L-iA/L^e-^ (4) 

Such a scaling for clusters up to L = 28 with one 
spin 3/2 in the system allows the determination of 
Aoo({X^i = 1/2, Ti = 3/2}) which remains positive up 
to very large values of J' (at least 2 J). Let us turn 
to the case of the first excited state of H{{Ti = 1/2}). 
Again, a scaling of the type = Aoo + {A/ L)e~ L ^ is 
expected with £ diverging and A^dTi = 1/2}) = at 
the isotropic point J' = J. Since the systems contains 
only spin 1/2 is that sector we performed diagonaliza- 
tions up to L = 48 to calculate A M . The results (see Fig. 
[|) clearly agree with what was expected even if the pre- 
cise behavior of A M ~ S 2 / 3 |L2] (up to logarithmic cor- 
rections) as the dimerization S = (1 — J'/J)/(l + J'/J) 
goes to zero is difficult to extract from numerical simu- 
lations. 

Discussion. The results are quite clear: For large sys- 
tems, J' c is larger than J, namely J' c ~ 1.15 J. Then to 
understand the low energy properties of the model in the 
parameter range J' < J and in the thermodynamic limit, 
one can use 7i({T, = 1/2}) as an effective Hamiltonian. 
This means that the low energy physics can be described 
by the dimerized, spin 1/2 Heisenberg chain with alter- 
nating exchange integrals J and J'. This model has been 
extensively studied in the context of spin-Peierls systems 
P^J . The main property, already used to perform the 
finite-scaling, is that there is a singlet-triplet gap in the 
system as long as J' < J, and that this energy gap closes 
at the point J' = J. Besides, the elementary excitations 
are boundstates of spinons due to the confinement poten- 
tial introduced by the dimerization |l3"| . Regarding the 
properties of the original problem, one should not for- 
get that all these states are degenerate, and in particular 
that there is, as in the previous case, a residual entropy 
per spin equal to 2 1 / 4 . 
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III. DISCUSSION 

The main motivation in undertaking the present study 
was to shed some light on the properties of the pyrochlore 
antiferromagnet. This system is a three dimensional 
structure that can be thought of as an array of coupled 
tetrahedra. The current situation as far as a theoretical 
understanding of that model is concerned is the following: 
There is a clear evidence that the spin-spin correlation 
functions are extremely short ranged, and there is some 
numerical evidence that there is a singlet-triplet gap in 
the spectrum . Both properties are indeed satisfied by 
the present ID model of coupled tetrahedra in the param- 
eter range J' < J. In particular, since the ground states 
can be written as products of local singlets, the spin- 
spin correlation functions indeed decay very fast. The 
fact that the gap closes when J' — J is not really in- 
consistent with the case of the pyrochlore since there is 
an important difference in the way tetrahedra are cou- 
pled in both cases: Pairs of tetrahedra are never coupled 
through more than one spin on each of them in the py- 
rochlore structure whereas one of the tetrahedra is cou- 
pled through 3 spins in the present model. So J' = ,/ 
in the present case somewhat corresponds to a stronger 
coupling than for the pyrochlore. 

More importantly, the presence of a residual entropy 
per spin in the model of the present paper suggests that, 
if there is a singlet-triplet gap in the spectrum, there 
should indeed be low-lying singlets within this gap. This 
should be particularly easy to detect if one goes away 
from the standard pyrochlore and considers a dimerized 
version of the same model with weakly coupled tetrahe- 
dra since the gap becomes larger in that limit while at 
the same time the splitting between the singlet states 
decreases. Work is in progress along these lines. 

To summarize, we have proposed and solved a one di- 
mensional analog of the pyrochlore antiferromagnet and 
proved that it exhibits at least some of the exotic physics 
one can hope to find in very frustrated magnets, namely: 
A singlet-triplet gap, and a lot of low-lying singlets. 
Given the relative simplicity of the model, it is the au- 
thors hope that some compound can be synthesized with 
this kind of structure. In any case, the very simple pic- 
ture of quite unusual properties provided by this model 
is expected to serve as a useful guide in the search of new 
experimental realizations of very frustrated magnets. 

The numerical simulations were performed on the Cray 
supercomputers of the IDRIS (Orsay, France). 
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